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Abstract: Coherent population trapping at zero magnetic field was 
observed for nitrogen-vacancy centers in diamond under optical excitation. 
This was measured as a reduction in photoluminescence when the detuning 
between two excitation lasers matched the 2.88 GHz crystal-field splitting 
of the color center ground states. This behavior is highly sensitive to strain, 
which modifies the excited states, and was unexpected following recent 
experiments demonstrating optical readout of single nitrogen-vacancy 
electron spins based on cycling transitions. These results demonstrate for 
the first time that three-level Lambda configurations suitable for proposed 
quantum information applications can be realized simultaneously for all 
four orientations of nitrogen-vacancy centers at zero magnetic field. 
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Impurity spins in solids are appealing as physical systems for quantum information process- 
ing, combining long decoherence times with the possibility of large-scale integration based on 
semiconductor processing technology [1]. Of the various impurities offering access to individ- 
ual electron spins, one of the most promising is the nitrogen-vacancy (N-V) center in diamond, 
which consists of a substitutional nitrogen atom next to a carbon vacancy. Electron spin co- 
herence times up to 58 /is have been measured at room temperature [2], and optical readout of 
the electronic spin state of a single N-V center has been demonstrated [3]. Controlled coupling 
between electronic and nuclear spins has been reported [4, 5], offering potentially much longer 
storage times, and controlled coupling with nearby nitrogen impurities is also possible [6]. 
Single -photon generation at room temperature [7, 8] has established the potential of N-V cen- 



ters for quantum communication. Based on such experimental results, schemes for quantum 
memories and repeaters [9] and quantum computation using electron spins [10, 11, 12] have 
been proposed. 

Most experiments conducted on N-V centers to date have relied on microwave fields to ma- 
nipulate the electron spins and have used optical fields mainly for readout. All-optical con- 
trol would allow for spatially selective addressing of single or a few N-V centers [10] and 
could be used in nonlinear-optical devices based on electromagnetically-induced transparency 
(EIT) [13, 14]. It is generally believed that the negatively charged N-V center has optical tran- 
sitions between the 3 A ground states and 3 E excited states at zero magnetic field that are almost 
perfectly spin-preserving. This situation would be advantageous for non-destructive readout of 
electron spins through photoluminescence detection, but — if this is the general case — poses a 
problem for all-optical control of the electron spins. All-optical control requires a A system, 
which has two ground-state levels connected to a common excited level by optical transitions. 
One solution is to work close to an avoided crossing of the m s = and m s = — 1 ground states 
that occurs at a particular magnetic field, as described in Ref. [15], where EIT was first reported 
in this system. However, previous spectral hole-burning studies [16, 17] suggest that the optical 
transitions are not always spin-preserving, even at zero magnetic field. In these experiments, 
anti-hole features (increases in photoluminescence) were observed when two excitation lasers 
were separated by ±2.88 GHz, the crystal-field splitting energy between the m s = and the 
degenerate m s = ±1 ground states. These features can only appear if one or more excited states 
couples to both the m s = and m s = ±1 levels. Thus there is an apparent contradiction between 
the appearance of these anti-hole features and the more recent experiments demonstrating op- 
tical readout based on cycling transitions. A possible clue in Ref. [17] is that the anti-hole 
structure showed a dependence on the center frequency, which could be interpreted as a strain 
effect. Nevertheless, the excited-state structure remains poorly understood, and it is important 
therefore to clarify the conditions necessary to obtain spin-nonpreserving transitions as needed 
for devices based on EIT and all-optical control of electron spins. 

Here, we report spectral hole-burning experiments on a sample with little inhomogeneous 
broadening and a range of strain conditions obtained by ion implantation. We find that it is 
possible to realize a A system even at zero magnetic field as confirmed through the observation 
of coherent population trapping [18, 19], the same mechanism responsible for EIT. The results 
confirm that strain has a strong effect on the excited-state fine structure [17] and suggest that 
N-V centers can be engineered to have either spin-preserving or spin-nonpreserving transitions 
as needed for a specific application. 

The sample used for this study was a type lb high-pressure, high-temperature (HPHT)-grown 
diamond crystal obtained originally from Sumitomo. Figure 1(a) shows a composite photolu- 
minescence image obtained under 532 nm illumination. The bright squares were implanted with 
2MeV He + ions at a concentration of ~ 10 16 cm~ 2 . The ion implantation breaks bonds [20]. 
Above a certain critical damage level, approximately 10 22 cm~ 3 vacancy concentration [21, 22], 
upon annealing the damaged material will relax to sp 2 bonded carbon. This graphitized material 
has a lower density than diamond which causes volume expansion and in turn produces strain 
in the surrounding regions. Annealing was performed at a temperature of 1400 °C for 15 min- 
utes in vacuum. Although nitrogen occurs throughout the crystal, most of the luminescence was 
observed to originate from close to the surface. Another important feature of this sample is the 
presence of distinct sectors. Growth sectors are regions that have grown from particular faces of 
the original seed crystal. Because different faces grow at different rates and incorporate impu- 
rities at different levels, the sample shows marked discontinuities at sector boundaries. In two 
particular growth sectors, one of which is outlined in Fig. 1(a), the inhomogeneous linewidth 
of the zero-phonon line (ZPL) for negatively charged N-V centers at 637 nm was approxi- 
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Fig. 1. (a) Composite photoluminescence (PL) image of a portion of the sample. The 
narrow-linewidth sector described in the text is outlined in green, (b) PL spectrum showing 
the zero-phonon line for negatively charged N-V centers, measured at location 1 in the im- 
age, (c) Polarized PL spectra from location 2. "H" and "V" are for excitation and collection 
polarizations oriented horizontally and vertically, respectively, compared with the image, 
(d) Schematic energy level diagram. 

mately 10 — 20 GHz, exceptionally narrow compared with the more typically observed value of 
~ 750 GHz [23]. Figs. l(b,c) show low-temperature photoluminescence spectra from two lo- 
cations, measured at 10K using 532nm excitation. At location 1 a single unpolarized peak was 
observed, while at location 2 the ZPL was split into two orthogonally polarized components. 
We attribute the splitting at location 2 to strain induced by the nearby implanted region, the 
edge of which is approximately parallel to one of the cubic crystal axes. By comparison with 
the stress measurements in Ref. [24], we estimate that the 44.5 GHz splitting corresponds to an 
equivalent uniaxial stress of 43 MPa. N-V centers can be oriented along four possible directions 
([111], etc.), but the splitting should be the same in each case for strain along [100]. A schematic 
energy-level diagram, based on Refs. [25, 26], is shown in Fig. 1(d). 

In the two-laser experiments, the zero-phonon line was excited on resonance using two 
continuous-wave, tunable external-cavity diode lasers. One of the lasers (laser 2) was held fixed 
in frequency while the other (laser 1) was scanned over a range of up to 80 GHz. The scanning 
rate was calibrated by means of a Fabry-Perot cavity. The laser stability was approximately 
1 MHz over 1 s and 10 MHz over a period of 10 minutes. A weak repump laser operating at 
532 nm was also applied continuously to re-activate N-V centers that bleach after many excita- 
tion cycles, an effect that is due perhaps to photo-ionization [27]. The lasers were focused by a 
0.6 numerical-aperture microscope objective through the cryostat window onto the sample sur- 
face with a spot size of approximately 3 — 4/im. The resulting photoluminescence (PL) was col- 
lected by the same objective, filtered to remove laser light, and sent to an avalanche-photodiode 
photon counter (Perkin-Elmer SPCM). A bandpass filter selected a 40nm bandwidth centered 
at 700 nm for detection of emission into the broad phonon sidebands. This provided a measure 
of the excited-state population. For most of the data shown below, laser 1 was scanned at a 
repetition rate of 1 Hz, and the results from many scans were summed. 

Results from typical two-laser measurements performed at low excitation power (0.5 /xW) 
at location 2 are shown in Fig. 2. With one laser fixed on resonance with the long-wavelength 
component of the ZPL, the frequency of the second laser was scanned across the entire ZPL 
structure. In Fig. 2a, a pattern of peaks appears prominently around the frequency of the fixed 
laser, and these peaks are shown with higher resolution in Fig. 2b. The pattern of peaks is 
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Fig. 2. (a) Two-laser scan at location 2, with the fixed laser tuned to the center of the long- 
wavelength peak in Fig. lc, producing the prominent anti-hole features. Both lasers were 
polarized along "V." (b) Higher-resolution scan of the anti-hole features, (c) Energy level 
diagrams explaining the individual peaks. The orange (red) arrows represent transitions 
driven by the fixed (scanning) laser. 



associated with the energy level structure and transition strengths at this particular location. As 
explained in Refs. [16, 17, 25], photoluminescence peaks occur at two-laser detunings for which 
each laser is resonant with a transition involving a different ground state. In these situations, the 
N-V center is excited from either ground state, producing continual photoluminescence. If only 
one ground state is excited, the population is driven to the other ground state (optical pumping), 
and photoluminescence is suppressed. 

The observed two-laser spectrum can be explained by an energy-level diagram with three 
excited states, shown in Fig. 2c, where the energy spacing between the ground and excited 
states is allowed to vary due to inhomogeneous broadening, but the level spacings within the 
ground and excited-state manifolds are constant. It is known that the ground states, formed 
from an orbital singlet, are relatively insensitive to stain [28], while the excited states, formed 
from an orbital doublet, are highly sensitive to strain [17]. In the excited states, strain lowers the 
symmetry, causing first a splitting of the orbital doublet [24], and then a modification of the spin 
sublevels through spin-orbit coupling [25, 26]. The peak positions in Fig. 2b thus indicate the 
excited-state fine structure. We observed that the peak positions and intensities in the two-laser 
spectra depend strongly on position relative to the edge of the heavily implanted region, from 
which we can infer a strain dependence. However, the peaks at detunings of ±2.88 ±0.01 GHz, 
matching the crystal-field splitting of the ground states, do not shift measurably with strain, 
and only their height changes. At these detunings the two lasers couple two ground states to 
the same excited state, and the presence of these peaks suggests that a A system with allowed 
spin-nonpreserving transitions is present. 

Observation of coherent population trapping requires only increasing the excitation power. 
Figure 3 shows two-laser spectra measured with 5 jJ.W excitation power at location 1 (Fig. 3a) 
and location 2 (Fig. 3b, c). The single-laser spectra in Fig. 3(d,e) indicate the corresponding 
spectral positions within the inhomogeneously broadened ZPLs. The narrow dips appearing 
within the 2.88 GHz anti-holes are the signature of coherent population trapping. Close to two- 
photon resonance, a "dark state" can form, which is a coherent superposition of ground states 
that has no net transition matrix element into an excited state due to destructive quantum in- 
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Fig. 3. (a,b,c) Two-laser scans showing anti-hole features and coherent population trapping 
effect at sample location 1 (a) and at location 2, measuring the long-wavelength (b) and 
short-wavelength (c) components of the zero-phonon line. Insets: higher-resolution scans 
of the —2. 88 GHz features. Excitation powers: 5 flW for the 637 nm lasers and 2/i W for the 
repump. (d,e) Single-laser scans at locations 1 and 2, respectively. Arrows indicate spectral 
positions studied in (a,b,c). 



terference [18, 19]. This feature only appears when the transition Rabi frequencies are large 
compared with the relevant decoherence rates. Observation of this feature demonstrates unam- 
biguously that we can find A systems with allowed spin-nonpreserving transitions. At location 
1, where no strain splitting is resolved, the effect was only observed on the long-wavelength 
side of the ZPL. On the short wavelength side, the anti-hole pattern was different, and the 
2.88 GHz anti-holes were almost absent. The random fluctuations in strain or electric field re- 
sponsible for the inhomogeneous broadening apparently has a strong effect on the excited-state 
fine structure and optical transitions [17]. The only distinct anti-hole at this location is the 
2.88 GHz feature associated with the ground-state splitting, suggesting a random variation in 
the excited-state structure. At location 2, where the anti-hole structure is clearly resolved, the 
2.88 GHz anti-holes are much more prominent in the long-wavelength ZPL component, and the 
coherent population trapping feature was easily observed there. For the short-wavelength ZPL 
component, these anti-hole features are much weaker, suggesting that the optical transitions are 
nearly spin-preserving, but the coherent population trapping feature could still be observed. 

At zero magnetic field it is possible for all four orientations of N-V centers to show co- 
herent population trapping behavior. To determine the orientation dependence, we performed 
two-laser measurements with a weak magnetic field (~ 100 Gauss) applied to lift the degener- 
acy of the m s = ±1 ground states. The magnetic field direction was adjusted so that each of the 
four orientations had a different energy splitting. As shown in Fig. 4, the photoluminescence 
dip on two-photon resonance splits into a maximum of eight components corresponding to the 
four orientations. For the low-strain region (location 1), the dips are divided into two groups 
responding to polarizations along [110] and [110]. We expect that one polarization primarily 
excites N-V centers oriented along [111] and [111], while the other polarization primarily ex- 
cites centers oriented along [111] and [111]. For polarization along [010], all four orientations 
produce photoluminescence dips. For the long-wavelength peak at location 2, eight dips appear 
for [010] polarization, showing that all four orientations can participate even in the strained 
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Fig. 4. Two-laser measurements with a weak magnetic field applied to lift the m = ±1 
degeneracy, (a) Measured at location 1 . Half-wave-plate (HWP) angle 40° corresponds to 
polarization approximately along [010]. (b) Measured at location 2 for polarization along 
[010] (V). All plots are scaled and shifted for clarity. The measured intensity decreases by 
approximately 8% for the dips in (a) and 4% for the dips in (b). 

case. The possibility of obtaining equivalent A systems for all four orientations is an advantage 
over the anti-crossing scheme in Ref. [15], where only one orientation participated. 

Finally, we show that the shape and excitation power dependence of the 2.88 GHz peaks can 
be explained by a theoretical model that takes into account the inhomogeneous broadening of 
the ground-to-excited-state transitions. The power dependence measured at location 1 is shown 
in Fig. 5. At low power, a shallow dip approximately 10MHz wide appears which then deepens 
and broadens with increasing power. To fit these data, we use a simplified model with three 
levels in a A configuration. The two ground states represent the m = (| 1)) and a superposition 
of the m = ±1 (|2)) levels of the ground state manifold, with a fixed frequency difference of 
2.88GHz. The excited state (|3)) is given a variable frequency following a normal (gaussian) 
distribution ((7 = 10 GHz). The model includes all possible population relaxation and dephasing 
terms under the simplifying assumption that the relaxation rates and transition strengths are the 
same for the two lower levels. We first find the steady-state density matrix for a single center and 
then average over the inhomogeneous distribution. The fitted curves, shown in Fig. 5, are pro- 
portional to the steady-state population of the excited state, plus a linearly sloped background 
added as an additional fitting parameter representing contributions from other subsets of levels. 
The theoretical curves match the data quite well, confirming that the observed behavior can be 
explained in terms of an inhomogeneous distribution of three-level systems. From the fits we 
estimate that the pure dephasing rate between the excited and ground states varies in the range 
713/27T = 2—11 MHz over the range of excitation powers used, and this should be compared 
with the radiative linewidth of 13. 4 MHz used in the fits. Most importantly, the fits can deter- 
mine the value of the ground-state decoherence parameter 712, since this parameter is closely 
related to the shape of the dip within the anti-hole feature. We find yn/2K = 3.5 — 7 MHz. 
Laser instability can contribute to the decoherence, but initial measurements using an electro- 
optic modulator in place of a second laser suggest that the true Y\2 is not much smaller than 
these values for this sample. Possible sources of ground-state broadening include inhomoge- 
neous variation of the ground-state crystal-field splitting [28], stray magnetic fields which can 
affect the four orientations differently, interactions with N impurities [6] and hyperfine coupling 
with 14 N and 13 C nuclei. The much longer coherence lifetime reported in Ref. [2] was obtained 
using spin-echo techniques to remove the effect of inhomogeneous broadening. 

The coherent population trapping effects reported here suggest that EIT and all-optical spin 
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Fig. 5. Coherent population trapping at various excitation powers at location 1 (same as 
in Fig. 3a). blue:data, red:fit. The indicated powers correspond to both 637 nm lasers. 
For excitation powers 0.5, 1.6, 5, and 16/1 W the fits used Rabi frequencies (same for 
both transitions): 10.1, 14.2, 17.9, and 22.0MHz; effective population decay rates between 
ground states: 0.1 1, 0.20, 0.35, and 0.83MHz; excited-state decoherence: 2.4, 3.0, 6.3, and 
11.4MHz; ground-state decoherence: 7.0, 5.8, 3.6, and 3.8MHz. 

manipulation should be possible with N-V centers even at zero magnetic field. A small amount 
of strain, which can be introduced through ion implantation or other fabrication methods [29], 
can determine whether spin-nonpreserving optical transitions are allowed. For strain along 
[100], all four orientations can produce a similar A system. Therefore N-V centers have a 
flexibility in their excited-state level structure that makes them suitable either for single-spin 
readout through photoluminescence detection or for optical devices based on EIT and Raman 
transitions. 
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